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Abstract 

Constrained orthogonal polynomials have been recently introduced in the study of the Hohenberg- 
Kohn functional to provide basis functions satisfying particle number conservation for an expansion 
of the particle density. More generally, we define block orthogonal (BO) polynomials which are 
orthogonal, with respect to a first Euclidean inner product, to a given i-dimensional subspace Si of 
polynomials associated with the constraints. In addition, they are mutually orthogonal with respect 
to a second Euclidean inner product. We recast the determination of these polynomials into a general 
problem of finding particular orthogonal bases in an Euclidean vector space endowed with distinct 
inner products. An explicit two step Gram-Schmidt orthogonalization (G-SO) process to determine 
these bases is given. By definition, the standard block orthogonal (SBO) polynomials are associated 
with a choice of £i equal to the subspace of polynomials of degree less than i. We investigate their 
properties, emphasizing similarities to and differences from the standard orthogonal polynomials. 
Applications to classical orthogonal polynomials will be given in forthcoming papers. 

PACS numbers: 02.10.Ud, 02.30.Gp, 02.30.Mv, 21.60.-n, 31.15.Ew, 71.15.Mb 
1. Introduction 

Recently, B Giraud et al. [1-3] have considered new sets of constrained orthogonal polynomials. Basically, 
these real polynomials Pn{x), n = 1,2, ... of exact degree n (i.e. the coefficient of x" is nonzero), satisfy the 
constraint of vanishing average ^ with a non-negative weight function w{x) on a real interval [a, b], 

b 

dxw{x)Pn{x) = n^l,2,.... (1.1) 
In addition, they are orthogonal on the same interval with a distinct non-negative weight function W2{x), 

b 

dxw2{x)Pm{x)Pnix) (X Sm.n m,n = 1,2,.... (1.2) 

The constant polynomial Pq is excluded since it does not fit the constraint. As a result, these polynomials 
does not form a complete set, but span a subspace which can be well-suited to specific applications. Thus, 
in various problems of mathematical physics, one considers an unknown function f{x) which has to satisfy 
a similar constraint, dx w{x) f {x) — 0. One way to take into account this constraint readily, is to expand 
f{x) in terms of the basis polynomials Pn{x), n = 1, 2, . . ., as f{x) — fnPn{x), with Uc a possible cut- 

off. The choice of the weight functions is dictated by the problem investigated and physical considerations, 
giving, possibly, a physical meaning both to the polynomials and to the possible truncation of the functional 
space. 

^ A constraint of 'vanishing momentum', dx w(x)xPn(x) = 0, is also considered in [1] section 3. See 
footnote 44. 



This method has boon apphcd to the Hohenberg-Kohn variational principle [4] for the ground stgroiind 
stateate energy, originally established for an interacting electron gas in a local spin-independent external 
potential v{r) leading to a non-degenerate ground state. The variable function is the particle density n(r). 
The Hamiltonian reads Ti. = T + V + U where T is the kinetic energy, V corresponds to the local one-particle 
external potential v{r) and U is the two-particle interaction (e.g., the Coulomb interaction). With number of 
particles N, and their mutual interaction U specified, let V be a set of local one-particle potentials ?;(r) such 
that, for each w(r), there exists a non-degenerate A^-particle ground state I^E") with energy E, that is solution 
of the Schrodinger equation = E]"^). From l^*), one calculates the ground state particle density n(r), 

see, e.g., equation (1.11) below. By definition, the number of particles is a functional of n(r) such that 



Thereby, one defines the mappings v{r) I^P) ^ n{r). Let the set J\f of particle density functions be defined 
as the image of V in the resulting mapping v{r) i-^ n{r). Since a real additive constant in v{r) has no 
effect on the ground state I'if) and thus on n(r), potentials in V differing by a real constant are considered 
equivalent. Let {w(r)} denotes the equivalence class. Then, the Hohenberg-Kohn theorem states that 
(i) for all v{r) in V and all n(r) in Af, there are one-to-one correspondences {^^(r)} \^) <-*■ n(r); 
(a) according to (i), let I^E") be the N -particle ground state corresponding to any n(r) in Af, then, 



is a universal functional, i.e. it does not depend on the external potential v{r) and is valid for any number 
of particles N = N[n]. Now, for any Vo{r) in V, the energy functional. 



assumes its minimum value on J\f for the exact ground state density no(r) corresponding to wo(r) according 
to (i) ^. Then, £'^(,[no] is equai to the corresponding N-particle ground state energy Eq = min„g_/v [n] . 
Let us note that the implementation of this variational principle raises several fundamental problems about: 
(i) the possible extension of TV to wider spaces of trial particle density; (ii) the existence of a functional 
derivative to formulate the variational principle as, S£vo[n]/Sn{r) =0 [5] ^. 

Based on the work of Hohenberg and Kohn and several extensions [5] ^, the density functional theory 
has become a standard approach to investigate the properties of quantum interacting many-particle systems 
in terms of the particle density n(r). This method has been frequently used in several branches of chemistry 
and physics, e.g., in atomic, molecular and nuclear theoretical physics. Along the Thomas-Fermi approach 
and the work of Kohn and Sham [6] (replacing direct variations with respect to the particle density by an 
auxiliary orbital picture), several functionals have been successfully used [5, 7]. Then, a constructive study of 
the ground state is provided by standard perturbation theories (e.g., particle- hole excitations, configuration 
mixing), around a mean field first order approximation. The success of these methods is based on the 
existence of suitable truncations to relevant subspaces. 

In view of their use in some problems of nuclear physics, B Giraud et al aim to make a similar approach 
in the space of particle density functions instead of the space of wavefunctions. In other words, can the 

^ The proof relies on the Rayleigh-Ritz minimum principle for the energy functional E[^] := in 
which the wavefunction 4" is the variable function. 

^ See, e.g., [5] sections 2.1 and 2.3 and the references therein, especially the works of Lieb and Levy. The 
particle density n(r) G J\f considered in the Hohenberg-Kohn theorem is said to be a pure state v-representable 
since it is the density of a pure ground state l^*} for a specific local external potential v{r) £ V. Now, for any 
given non-negative n(r), normalized according to equation (1.3), does there exists a local external potential 
v{r) so that v{r) n(r)? The answer can be negative, even when the formalism is extended to non-pure 
N-particle states described by density operators. 

^ Degenerate ground state, finite temperature ensembles, etc.; see, e.g., [5] section 3 and the references 
listed in this review book. 




(1.3) 



Vn(r) e AT, n(r) ^ \^) 



F[n\ := {^\T + U\^) 



(1.4) 




(1.5) 
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functional spac;c of n(r) be trTincatcd to a subset of meaningful 'density modes'? [1 3]. An ultimate 
goal would be a constructive approach to the Hohenberg-Kohn functional. As a first step, the one-to-one 
correspondence between the particle density and the external potential (up to within an additive constant) 
is investigated through the related values of Sn{r) and Sv(r) around the true solution no(r) for a given 
Vo{r). Then, the question is: can the functional spaces for 5n{r) and for Sv{r) be defined by 'suitable' 
basis functions, such that the truncation of theses functional spaces to trial subspaces of few basis functions 
(or 'modes') be relevant? In heuristic approach, forgetting the questions brought up previously about the 
functional spaces, we focus on the following basic constraints: (i) n(r) has to be non-negative and the particle 
number conservation (1.3) yields 

'rfr^n(r)=0; (1.6) 



(ii) the one-to-one mapping is between {w(r)} and n(r), therefore 6v{r) must not be a nonzero constant, 

Sv{r) ^ const . (1- 7) 

A standard way to take into account the particle number conservation with the variational formulation is 
to introduce a Lagrange multiplier /x such that (5(iSt,„[n] — /U / dr n{r)) / 5n{r) = 0. Once again, this requires 
an extension of the functional space to particle density normalized to a not-necessarily-integer value as in 
equation (1.3) ^. The method sketched in our first paragraph presents an alternative overcoming these 
difficulties by considering trial variations Sn(r) which inherently satisfy the constraint (1.6). This is done 
as follows in a simple toy model. As a usual first approximation in the nuclear shell model, a system of N 
independent fermions in a one-dimensional harmonic oscillator potential vo{x) := ^ x"^ is considered in [1-3] 
^, neglecting the two-body interactions. Then, the single-particle orbital wavefunctions are defined by [8] '', 



1/^ 

2V (ix2 



x^ ] h -ipj (x) = ej ipj (x) (1.6 



i;j{x) ={2"n\)-'^w-ie-i-'Hjix) = ^ + j j = 0,1, . . . (1.9) 

where Hj{x) is the Hermite polynomial. The TV-particle ground state wavefunction 5'o(a;i, . . . ,xn) is the 
properly normalized N x N Slater determinant built from the N lowest energy single-particle eigenf unctions, 

*o(a;i,...,a;jv) = {Niyi det[i)j{xk)] j=o.....n-i . (1.10) 

Je = l,...,iV 

The particle density, satisfying the normalization (1.3), reads 

/oo /"OO N—1 

dx2'-' / fix AT 1^0 (^,^2, ■ ■ •,xn)\^ = ^ i^A^Y . (1.11) 



3=0 



Now, it follows from the first order perturbation theory [8] ^that the response to a variation Sv{x) is for 
each single-particle wavefunction, 

dijjix) = Vj(a;)-^^^^^^ ii^jM^Jj) = r dx^j{x)5v{x)ijj{x) (1.12) 



J=N 



where j = 0, . . . , TV — 1 and J = N, ... ,00 are the hole and particle indices, respectively. As a result, from 
equation (1.11), the variation of the particle density is ^ 



N-l 00 

Sn{x) = 2 V V ,Pj{x)Mx) (1.13) 



j=0 J=N 



5 



7 



See, e.g., [5] section 2.4 and [7] section 1.4.5, the problem of derivative discontinuities. 
See [1] section 4, [2] section 4 and [3] section 2. 
See, e.g., [8] appendix B (B.65) and (B.70). 
8 See, e.g., [8] chapter XVI, section 2 (XVI.16). 

® The positivity of the density po + Sp is guaranteed at first order, since dp follows from a first order 
calculation of the perturbed N-particle wavefunction -|- 5^. 
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thereby, providing the functional mapping Sv(x) Sn{x), hnear in first order perturbation theory. As usual, 
for practical study, this functional correspondence can be transformed into a discrete (possibly infinite) 
linear problem by expanding both 5n{x) = Y^^=o^''^jVj{^) 'i^id 5v{x) = Y^^=o^''^k4'k{x) on suitable bases, 
{iPn{x), n = 0, 1, . . .} and {(j)n{x), n = 0, 1, . . .}, respectively. Now, since tpj{x)'ipj{x) oc exp(— x^) times a 
polynomial, Sp{x) can be expanded readily on functions ipn{x) ■= exp(— a;^)P„(a;), with {P„(a;), n = 0, 1, . . .} 
any basis of real polynomials. It is convenient to choose these polynomials of degree n and such that the 
functions i^„(a;), n = 0, 1, . . . be orthogonal, 

dx i^rn{x)tpn{x) = / dx t'"^"^ Pm{x)Pn{x) <X. 5m,n ■ (1-14) 
-oo J — oo 

This amounts to choose the polynomials Pn{x), n = 0, 1, . . . to be orthogonal with respect to the weight 
function W2{x) := cxp(— 23;^). In addition, the particle number conservation (1.6) will be satisfied trivially 
if each basis function (fin{x) fulfils the constraint, 

/oo fOO fOO 

dxtpn{x)= / dx e~''^ P„{x) (X / dx e'''^ PoPnix) = (1.15) 
-oo J — oo J — oo 

i.e., if the polynomial Pn{x) is orthogonal to the constant polynomial Po with the weight fimction w{x) := 
exp(— x^). As already noted, this excludes the value n = 0. Thus, for this toy model within the framework 
of perturbation theory and particle-hole excitations, the functions (pn{x) = w{x)Pn{x), with polynomials 
Pn{x), n = 1, 2, . . . satisfying equations (1.1) and (1-2), assuming such polynomials exist, provide a functional 
space for 5n{x) taking into account the particle number conservation, essential for the Hohenberg-Kohn 
variational principle. Now, what about the choice of the basis functions (t)n{x) for the expansion of 5v{x)l 
Although the constraints (1.6) and (1.7) arc different, it is convenient to choose the same basis, i.e. (t>n{x) '■= 
ipn{x), n = 1,2,... (otherwise, connection coefficients between the two bases would be required in the 
calculation). Nevertheless, it should be noted that the constraint (1.6), J^^dxSv{x) = 0, implies Sv{x) 
satisfies the constraint (1.7), but this not a necessary condition and therefore, this corresponds to a 
particular restriction of the functional space for Sv{x). Then, it is argued in [1-3] that these basis functions 
are good candidates to define relevant particle density modes. The same approach can be generalized for any 
one-particle potential wo(r), not only in a one-dimensional model. Thus, in connection with the Laguerre 
polynomials, the following basis functions and weight functions are considered in [2] 

tfinix) := e-^P„(a;) w{x) := e'^'x'^-^ W2{x) := e'^'^x'^-^ . (1.16) 

The functions (Pn{x), n = 1,2,... have to satisfy the orthogonality relations (1.14) and (1.15) with now 

X G [0, oo) as, e.g., the radial variable in a d-dimensional space. 

The purpose of this paper is to make a systematic study of this kind of polynomials defined by equations 
(1.1) and (1.2). The problem can be generalized as follows. For given nonzero and non-negative weight 
functions w{x) and ^U2{x) on a given interval, and for given i linearly independent polynomials ei{x), . . . , ei(x) 
(possibly associated with i 'constraints' ^^), spanning a subspace £i of a space £ of polynomials: (i) does 
there exist a subspace £2 orthogonal to £1 with respect to the weight function w(a;), together with £1 and £2 
complementary in £7 (ii) is it possible to define, and compute, within £2 an orthogonal basis with respect 
to the other weight function W2(x)? Thereby, one defines what we call block orthogonal (BO) polynomials, 
instead of constrained orthogonal polynomials, to underlined the fact that the linear constraints are defined 

B Giraud et al argue in [1] section 1, [2] section 1 and [3] section 2, that 'the basis must be orthogonal 
to a flat potential' (i.e. constant). This is not true, e.g., 6v{x) = (poix) oc w{x) does not fit the constraint 
(1.6), but satisfies the constraint (1.7). A basis excluding a constant function could be {x", n = 1, 2, . . .}. 

For example, in atomic and molecular physics, one may think about the Coulomb interaction vo{r) oc 
1/x, where x is the radial variable in a 3-dimensional space. The single-particle orbital wavefunction includes 
an exponential term linear in x, see, e.g., [8] chapter XI, section 6. Then, it would be clumsy to consider a 
basis function ipn{x) with an exponential term quadratic in x, as for the harmonic oscillator. 

For equation (1.1), i := 1 and ei{x) := Pq = const. 
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by an inner product. Although tliv problcirn was motivated by polynomial considerations, it is worthwhile 
to point out what can be generalized and ascribed to Euclidean vector space and to polynomial algebra, 
respectively. 

This paper is organized as follows. The general problem in linear algebra is considered in section 2. 
For the sake of clarity, elementary results for Euclidean spaces are recalled in section 2.1 (see, e.g. [9-11]). 
Two BO subspaces and bases are defined and constructed in section 2.2. The well-defined method is based 
on a two step Gram-Schmidt orthogonalization (G-SO) process. The possible extension to more than two 
subspaces is discussed in section 2.3. These formal algebraic considerations are applied to polynomial vector 
spaces in section 3. After generalities about Euclidean vector spaces of polynomials in section 3.1 (see, e.g. 
[12 14]), BO polynomials for two subspaces arc constructed in section 3.2. Then, we focus, for a given 
non- negative integer i, on the real BO polynomials Pi-^n{x), n = i,i + 1, . . . of exact degree n, orthogonal to 
the polynomials of degree less than i with the weight function w{x), and mutually orthogonal with respect 
to the weight function W2{x) 



/ 

J a 



b 

dx w{x)x"'' Pi-nix) = m = 0,...,i-l n = i,i + l,... (1-17) 



/ dx W2{x)Pi.rnix)Pi-n{x) oc dm,n m,n = i,i + 1, . . . (1-18) 

J a 

generalizing equations (1.1) and (1.2) for i > 1. We call these polynomials standard block orthogonal (SBO) 
polynomials. Their general properties are investigated in section 3.3, comparing them with the properties 
of the standard orthogonal polynomials. We give our conclusions in section 4. For completeness, and to 
underline the similarities to and the differences from the study of BO polynomials, definition and properties 
of standard orthogonal polynomials are recalled in appendix A. Examples of BO polynomials for three 
subspaces are given in appendix B. 

Throughout the remaining of this paper, the following conventions and notations are used: 

- i, j, k, £, m, n and N denote non-negative integers; 

- a null sum is interpreted as zero while a null product or a null determinant is interpreted as unity; 

- matrices are denoted in boldface, e.g., A := {Aj^k)j,k=---', 

- for any function / : x f{x) (especially a polynomial), / stands for f{x); 

- monic polynomials, i.e. with the coefficient of the highest power equal to one, and also any related quantities 
are denoted by hatted letters, e.g., P„ and /i„; 

- in the differentiation and/or recurrence formulae, the polynomials with negative degree are set equal to 
zero, e.g., P_2 = P-i := 0; 

- standard orthogonal polynomials are denoted by Qn, n = 0, 1, . . ., where Q„ is of exact degree n; 

- classical orthogonal polynomials (e.g., the Hermite and Laguerre polynomials) are defined according to 
[14 16] 

- BO polynomials are denoted with a capital P, e.g., Pi-n for SBO polynomials. 
2. Block orthogonal subspaces and bases 

2.1. Basic definitions and properties of complementary and orthogonal subspaces 

Let £ be an iV-dimensional vector space over the real field R. A symmetric and positive-definite (thus 
non-degenerate) bilinear form on £, defines an Euclidean inner product or scalar product ( , ) such that 

yu,ve£ (u , v) = {v , u) e R (w,?i)>0 and (u , u) = <^ u = . (2.1) 

1^ In other words, the i linear constraints are associated with the i linearly independent polynomials 
ej-i{x) := x^ , j = 0, . . . ,i — 1, defining the subspace Vi :— £i of polynomials of degree less than i. 
" See, [14] chapter X, [15] chapter 22 or [16] section 8.9. 

1^ The bilinear form is said to be: (i) non- degenerate if for every u, {u , v) = Q implies v = 0] (ii) definite if 
(m , m) = implies u = 0; (iii) positive if {u , u) is positive for every nonzero u. See, e.g., [9] sections 59-61, 
[10] chapter XIV, section 1 or [11] 3.1-1. 
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Any A^i-dimcnsional subspace £i of £, not {0} or £ itself, has an infinity of cornplement,ary subspaces, all 
of dimension N2 = codim^ £1 := N — Ni, i.e. the codimension ^^of £1 with respect to £. Indeed, any basis 
{ei, . . . , ejvi} of £1 can be completed with vectors ei,. . . ,sn2 such that {ei, . . . , e^, ei, . . . , eatj} is a basis 
of £. Then, £1, . . . , 6^2 generate a complementary subspace £2 of £1. For any arbitrary vectors ui, . . . , u^^ 
in £1, not all zero, the N2 vectors, e'^ := Uj +Sj, j = 1, . . . , N2 generate another complementary subspace £2 
of £1 . Now, it is known that there is a unique subspace £2 ~ £1 such that, £1 and £^ are complementary 
and orthogonal with respect to the Euclidean inner product ( , ), 

£ = £ie £]r and {£1 , £^) = ^ Vu e 5i, V?; e {u,v) = 0. (2.2) 

Tien, £^, is called the orthogonal (with respect to ( , )), complement of £\ ^"^ . 

A constructive method to define £^ is to make use of the Gram-Schmidt orthogonalization (G-SO) 
This is a canonical inductive procedure of getting an orthogonal basis {Ei, . . . , En} of £, starting with an 
arbitrary basis {ei, . . . , ejv}- For bjj, j = 1, . . . , n some arbitrary real nonzero finite factors, one defines 
recurrently 

El := ei&rj 

E2 := (-£;i6i,2 + e2)b2^l 



Er, := (-£i6i,„ En-ibn-i,n + e„)6-'„ n = l,...,N (2.3) 

such that, 

n n 

{Ej , Ek) = hjSj^k j,k = l,...,n E„ = emam,n e„ = ^ Emhm,n an,n = Kjn ■ (2-4) 

m=l m=l 

Once a basis is chosen, the inner product is characterized by the metric tensor components, 

9j,k {cj , Ck) i,j = 1,...,N . (2.5) 
It is known ^°that En can be written formally in terms of a determinant as. 



7-^ ^n,n 1 . 

En = det 



(efc) 



n=l,...,N (2.6) 



3 = 1.- 
fe=l,. 



where Z„ is the nonzero Gram determinant 

Zo:=l Zn:=det[gj,k]j^k=i,...,n n = l,...,N. (2.7) 
See, e.g., [17] A, II, section 7, no 3, definition 2, p 99. 

The uniqueness of follows from the non-degenerate character of the inner product. See, e.g., [9] 
section 62 and the Projection theorem in section 66, [10] section 1 proposition 2 or [11] 3.3-3 and 3.3-4. In 
[8] chapter VII, section 4, complementary and orthogonal subspaces are said to be 'complementary'. 
18 See, e.g., [9] section 65, [10] chapter XIV, section 7, [11] 3.4-6 or [14] 10.1 (7). 

1^ If the inner product were not definite, an arbitrary basis vector as ei might be such that (ei , ei) 
vanishes. But the non-degenerate character of the inner product implies that any basis vector, say -Ei- of an 
orthogonal basis has to be such that {Ei , Ei) ^0, otherwise this nonzero vector would be orthogonal to all 
vectors in contradiction with the hypothesis. 

^'^ This is usually considered within the framework of orthogonal functions and especially polynomials. 
See, e.g., [12] (2.2.6) or [14] 10.1 (9), (10), 10.3 (3) and (4), with A;„ := 6-i„, G„ := Z„ and Cj+k ■= 9j,k- 

^1 The inner product being positive definite, the real symmetric matrix {gj,k)j,k=i,...,N is positive definite. 
All its principal minors are positive, see, e.g. [14] section 6.2.4. Hence Z„ 7^ for n = 0, . . . , -/V. 
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Actually, ail the relevant quantities can be expressed in terms of determinants with gj^k as entries, except 
possibly for the last row. One has for n = 1, . . . , A'', 



det 



{9n,k) 



3 = 1, 



n-1 



^n-1 



■det 



(.9j,fc) 



3 = 1,. ..,n- 
fc=l, . . . ,n 



= (-ir+"|^det[ff,,.].^....,„- 
z„_i 



(2.8) 
(2.9) 



It follows from equation (2.3) that for m = 1, . . . , n — 1, the orthogonality relation {Em , En) = yields 
hmbm.n — (Em 7 e„). Then, using equations (2.6) and (2.8) and the symmetry property, gk,n — 9n,ki 

the 

connection coefficient bm,n reads for 1 < m < n < N , 

1 bm,m / J , 
Om,n = T y (det 

This formula is not given in the standard textbooks already quoted. 

There is an infinity of orthogonal bases of the subspace generated by the n linearly independent vectors 
d, . . . , e„ The G-SO process shows that: for any set of nonzero finite factors {bjj 7^ 0, j = I, . . . , n}, 
there is a unique sequence of orthogonal basis vectors {Ei, . . . , £"„} satisfying equation (2.4), i.e. such that 
for m= 1, . . . , n, Em is a linear combination of ei, . . . , with the nonzero component on e^. 

Now, assuming {ei, . . . , cat^} is a basis of the subspace £1 and applying G-SO up to n = Ni, then 
{El, . . . ,E]\[^} is an orthogonal basis of £1. Carrying on this procedure up to n = N, the N2 orthogonal 
basis vectors E^^+i, . . . , En determine the unique subspace £1 satisfying equation (2.2). 

In the commonly used Dirac notation each vector t; in f is denoted as the ket \v). Each element 
of the dual space £* is denoted as a bra {u\. Then, the orthogonality and the closure relations for the basis 
{El, . . . , En} read 

N 

{Ej\Ek) = hjSj,k j,k = l,...,N Y.\^n)h-^ {En\ = h (2.11) 

n=l 

where Ig denotes the identity operator on £. The projector Xlf^ onto £1 along £1 , and the projector Ii^±. 

onto £1 along £1 are said to be orthogonal projectors or perpendicular projectors with respect to ( , ). One 
has 

JVi N 

ne, = Y.\En)h-^{En\ ns^= \En)h-'{En\ He, + 11^. = Is . (2.12) 

n=l n=Ni+l 

Although often convenient, this compact notation must be handled with care in the continuation, since several 
inner products are going to be considered simultaneously, see equation (2.20). Therefore, the Dirac notation 
will be used henceforth only when dealing with orthogonal projectors. 

Remark. Following the idea of 'constraint' considered by B Giraud et al for equation (1.1), or more generally 
(1.17), £2 = £1 can be alternately defined as the subset of vectors v m £ satisfying the i = Ni constraints, 
[cj , v) = 0, j = 1,. . . ,Ni, where ej, j = 1,. . . ,Ni are linearly independent vectors (a basis of £1 in our 
first formulation). Since these constraints are linear in v, £2 is a subspace. More precisely here, they are 
associated with the linear forms, {ej\, j = 1, . . . , A'l, and read {ej\v) = 0. It follows that v belongs to the 
kernels of each of these forms, each kernel being of codimension one Now, the non-degenerate character 

The sequence {E[, . . . , E'^} is also an orthogonal basis with the same normalization (2.4) if and only if 
it is obtained from {£^1, . . . , En} by an isometric linear mapping, E'^ := J2k=i Ekh^^^^Okjhy^ , j = 1, . . . ,n 
where, {Okj)k,j=i,...,n is any n x n orthogonal matrix in the orthogonal group 0(n). 
23 See, e.g., [8] chapter VII. 

2^ Due to the non-degenerate character of the inner product, the kernel of {ej\ is the {N — l)-dimensional 
hyperplane orthogonal to ej with respect to ( , ) in the A''-dimensional space £. 



[9j,k) 
(ek) 



■det 



[9j,k) 

{9n,k) 



(2.10) 
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of the inner product implies that these Ni forms are Unearly independent. Consequently, £2 is a subspace 
of codimension Ni . 

In a more general problem, independently of the existence of an inner product, £2 could be defined as 
the intersection of the kernels of Ni linearly independent linear forms fj,j = l,..., Ni in £* . For any given 
basis {cj, j = 1, . . . , N} of £, these forms are characterized by the coefficients fj^k '■= fji^k), j,k = 1, . . . ,N. 
Then, with v = J2j=i''^j^j' the subspace £2 is determined by the Ni linear equations with N unknowns, 
'^k=i fj,hVk = 0, j = 1, . . . ,Ni. The rank of this system being A''i, the solution v, defining £2 uniquely, 
depends linearly on A'^ — Ni parameters, e.g. wjvj+i, . . . , w^v assuming dct[fj,k]j.k=i....,Ni 7^ 

Nevertheless, the physical applications quoted in section 1 involve constraints which are particular linear 
forms, i.e. bras associated with an Euclidean inner product ( , ). Only this case is considered in this paper, 
studying below the interplay of the initial inner product (, ) with new inner products, (, )i on fi, and 
mainly ( , )2 on £2. 



2.2. Two block orthogonal subspaces and bases 

2.2.1. Definitions and properties. From above, for any subspace £1 of an Euclidean space £, with the inner 
product ( , ), the complementary and orthogonal subspace £1 is uniquely determined. Now, each of these 
subspaces can be endowed with a new Euclidean structure by any new Euclidean inner product ( , )i on f 1 
and ( , )2 on £^ (these inner products can be possibly defined on the whole space £). Then, £1 and £^ are 
said to be block orthogonal (BO) subspaces. 

Every vector in £ having a unique decomposition on the complementary subspaces £i and £1, these 
inner products induce a new Euclidean structure on £ with the inner product ( , )o defined by 

\/u = ui + U2&£ uie£i, U2(^£t I N / N ,/ s, /o1Q^ 
w , ^ c-L (m, -i^jo := (wi , t^iji + (W2 , ■^2)2 (2.13) 

or, equivalently, £\ and £^ being orthogonal with respect to ( , ), 

(u,v)\ Vu, u G £1 
(u, 'y)o := <( (u, w) = Vu £ £1 Vw G (2.14) 
[u , 11)2 Vw, 11 ^£^ . 

Now, an orthogonal basis {^i, . . . , (?ijv} of £ with respect to ( , )o, i.e. such that, 

{(t)j , A)o = HjSj,k j,k = l,...,N (2.15) 

or, from equation (2.14), 

, Ml = Hj5j^k 3,k = l,...,Ni 
'jAN^+k) =0 j = l,...,7Vi k = l,...,N-Ni (2.16) 

. (0jVi+j , 4>Ni+k)2 = H^i+jSj^k j,k = l,...,N -Ni 

can be obtained using G-SO again, separately in £1 and £1. Then, {^1, . . . , (I)Ni} and ■ ■ ■, <Pn} are 

said to be BO bases. 

It should be noted that 

- the inner product ( , )o defined on £ by equation (2.13) depends on £1 and on the three inner products 
(,),(, )i and (, )2; 

- an orthogonal basis of £1 with respect to ( , )i, {(pi, . . . , (pNi}, depends only on this inner product. It can 
be determined in the standard way, using G-SO to orthogonalize any basis {ei, . . . ,eNi} of £1 (possibly an 

When the linear forms are defined by /j := {ej\, j = l,...,Ni, as in the case we consider, then 
fj,k = 9j,k '■= (cj , Cfe) and thus Z^^ ^ 0, see footnote 21. 
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orthogonal basis with respect to ( , )). Thereby, with arbitrary real nonzero finite factors 6„,„, n= 1,. . . ,Ni, 
one has 

{ei,...,eArJ ' {(pi, . . . ,(j)Ni} 

9i;j,k ■■= {ej , ek)i {(t)j , (pk)i = Hj6j^k j,k=l,...,Ni (2.17) 

applying equations (2.4) (2.10) with 0„, _ff„ and giyj.k instead of /i„ and gj.k^ respectively; 
- on the other hand, an orthogonal basis of £^ with respect to (, )2, {(/'atj+i, . . . , <^jv}, depends on £i 
(regardless to the inner product (, )i), and on both inner products (, ) and (, )2- It can be obtained 
applying G-SO twice: 

(i) once to determine itself by the orthogonal basis with respect to ( , ), {Em^+i, ■ ■ ■ , E^}- With arbitrary 
real nonzero finite factors 6jVi+ri,iVi+n, n = 1, . . . ,N — Ni, one has readily from equations (2.4)-(2.10), 

{ei,...,e7v} — ^ > {El, . . . , E^n Emt^+1, ■ ■ ■ , En} 

gj.k ■■= (cj , ck) (Ej , Ek) = hj6j,k j,k = 1, . . . , N ; (2.18) 

(ii) once again to orthogonalize {i^jvi+i, • ■ • , E^} with respect to ( , )2- Thereby, with arbitrary real nonzero 
finite factors f3N-^+n,Ni+n, n = 1, . . . ,N — Ni, one has 

{ENi+i,---,E]sr} ^ > {4>Ni+i, ■ ■ ■ ,4>n} 

g2;j,k ■= {Eni+J , ENi+k)2 i4>Ni+j , 4'Ni + k)2 = HN^+jSj^k j , k = 1 , . . . , N - Ni (2.19) 

applying equations (2.4) (2.10) with (j)j<f-^^n, -ffwi+n and g2:j.k instead of £"„, /i„ and g^ fe, respectively. 

There is an inGnity of BO bases. They are related by isometric linear mappings within each one of the 
subspaces £i and S^-, see footnote 22. 

2.2.2. Projection operators. To use the Dirac notation, one has to take care of the dependence of the one-to- 
one mapping between the vector space and its dual space upon the inner product considered. Thus, to each 
ket \u), there corresponds now four distinct bras, denoted {u\, o{u\i i{u\ and 2{u\, and defined by, 

V\v)e£ {u\v) := {u , v) j{u\v) := {u, v)j j = 0,1,2. (2.20) 

The projectors orthogonal with respect to ( , ), Hg-^ onto £i and Ilg± onto £^ , can be readily expressed only 
in terms of bases of £ which can be split into a basis of £i and a basis of £i . This is the case for the basis 
{El, . . . , En} orthogonal with respect to ( , ). Then, the projectors are given by equation (2.12) in terms of 
the bras ( |. This is also the case with the BO basis {ipi, . . . , (f>N}. Then, with the bras o( | associated with 
( , )o, the orthogonality relation (2.15) and the closure relation read 

N 

o{ct>M = HjSj,k j,k = l,...,N ^|0,)iJ-io(.^,|=7f . (2.21) 
The projectors Hs^ and II^j. are given by 

Afi N 

n^, =^ !</.,) iIrio(</,^. I n^. = J2 \<t>j)Hi' o{<t>j\- (2.22) 

i=i j=Ni+i 

Now, one may ask to express these projectors in terms of the inner product ( , )2 and one of its orthogonal ba- 
sis. Assuming the inner product ( , )2 is defined on the whole space £, let us choose the basis {£^2;i, • • • , E2-n} 
one gets applying the G-SO procedure on the basis {E\, . . . , Ej^}, i.e. such that, 

G-SO 

{El,..., En} * {E2-1, . . . , E2-n} {E2;j , E2-k)2 = h2;jSj^k j,k = l,...,N 

n n 
E2;n = ^ Ema2-m,n En = ^ -E2;m&2;m,n n = 1, . . . , N (2.23) 
m=l m=l 



where, a2-n.n = b'zili,!!- '"' = I, ■ ■ ■ ,N are some arbitrary real nonzero finite factors. Then, the orthogonality 
and the closure relations in £ read 



N 



2{E2;j\E2;k) = h2;jSj,k j,k=l,...,N J2\^^'^^^2-^2{E2;n\=Is- (2.24) 

n=l 

Prom the closure relations (2.11) and (2.24), one has 

n 

\E2;n) = ^ \Em) (i?m|i?2;n) h^^^ {Em\E2;n} = Ct2;m,n (2.25) 

n 

\En) = ^ |-^2;m) /l2;TO 2{E2;m.\En) /12;to 2{E2;m\En) = &2;m,n (2.26) 



m=0 



where, using equation (2.23), the upper bound of the sums over m are set equal to n instead of A^. In other 
words, with the choice of basis {E2a, ■ ■ ■ , E2:n}, {Em\E2:n) and 2{E2;m\En) vanish for m > n. Finally, using 
equations (2.12) and (2.24)-(2.26), the projectors and II^j. can be expressed in terms of the bras 2( | as 



n N Ni 



= E E 1^2;,) b2;j,na2;n,k)h2.), 2{E2;k\ (2.27) 

j=l k=n 71=1 

N 

= E l-^^;™) 2{E2-n\ ' ^S, ■ (2.28) 

n=l 

2.3. Extension to more than two subspaces 

2.3.1. First generalization. One may ask wether it is possible to consider three, or more, BO subspaces. 
From above, for any given A^i-dimensional subspace £i of the Euclidean iVi.2-dimensional space £1,2 with 
the inner product ( , )i,2, there is a unique complementary and orthogonal (A^i,2 — iVi)-dimensional subspace 
£2 such that, according to equation (2.2), £1^2 = £\® £2 and {£\ , £2)1,2 = 0. Now, let us assume £1,2 is a 
subspace of an Euclidean TV-dimensional space £ with the inner product ( , ) . Iterating the procedure above, 
there is a unique complementary and orthogonal (iV — iVi^2)-dimensional subspace £^3, such that £ = fi,2©£3 
and {£\ 2 , £3) = 0. Thereby, both £1 and £2 are orthogonal to £s with respect to the same inner product 
(.)• 

More generally, for any given subspace £1 (B £2 of £ (e.g., such that {£1 , £2)1, 2 = 0) and two distinct 
Euclidean inner products ( , )i,3 and ( , )2,3 defined on £, does there exist a subset £3 such that, 

£ = £^(B£2®£3 and j.^f'f^''^^ (2-29) 

L (02 , ^3)2, 3 = \) c 

From the linearity of the inner product, if it exists, £3 is a subspace. Let us try to determine £3 by constructing 
one of its basis. Let {ei, . . . , ejv} be any basis of £ such that {ei, . . . , ejvi} and {ejvi+i, . • . , ejVi+jVa} are 
bases of £1 and £2, respectively. Now, setting 

N1+N2 

£„ := ^ emCm,n + eNi+N2+n n = I, . . . , N - {Ni + N2) (2.30) 
m=l 

these vectors are linearly independent. For n given, the Ni + N2 coefficients Cm,n, m = 1, . . . ,Ni + N2 have 
to satisfy the Ni + N2 orthogonality conditions of with £1 and £2 according to equation (2.29), 

f (ej , £„)i,3 = j = l,...,Ni ^^^^^ 

\ (CiVi+j , £n)2,3 = j = 1, . . . , 7V2 • 
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With the metric tensor components, ga,3;j,k '■= 



{ej , ek)a,3, j,k = l,...,N for a = l,2 and. 



\ fc=l,...,iVi+iV2 / ^ ' — i,---,JV2 / 

the Hnear equations (2.31) read in matrix notation, G (cm,„)^_^ N1+N2 ~ Now, det G is not a 

Gram determinant, as it would be if (, )i3 = (, )2,3. Therefore, G may be a singular matrix, and the 
linear equations (2.31) have a solution if and only if the rank of G and the rank of the so-called augmented 
matrix (GVn), obtained by adjoining to G the column Vn, are the same, say r. Then, the solutions 
for £„ depend linearly on Ni + N2 — r parameters. Since the principal minors det[gi^3-j^k]j,k=i,...,Ni and 
det[(j'2,3;7Vi+j,Afi+fc]j,fe=i,...,Af2 are Gram determinants, one has sup(Afi,7V2) < r < Ni + N2. Thereby, the 
problem of the existence of £3 satisfying equation (2.29), may have a unique solution, no solution or an 
infinite number of solutions. That these three cases do happen has still to be proven. In other words, 
it has to be checked that the relations required between the two metric tensors are compatible with their 
symmetric and positive-definite character. Examples with polynomials and Ni = N2 = I and N = 3 are 
given in appendix B. In the generic case, the rank of G is iVi + A^2 , i-c. det G 7^ 0, and £3 is defined uniquely 
by the basis vectors £1, . . . ,eN-Ni-N2- Thus, in this case, for given £1 C £1^2 C £ and the Euclidean inner 
products (, )i_2, (, )i,3 and (, )2,3, there are unique subspaces £2 and £3 satisfying equation (2.29) with 
£1,2 = £1® £2- Now, endowing £1, £2 and £3 with new Euclidean inner products ( , )i,i, ( , )2,2 and ( , )3,3, 
respectively, these three subspaces are said to be BO. A new Euclidean structure on £ is induced with the 
inner product defined as in equation (2.14) by, (u , v)o := (u , w)cr,cr if u, w G £0-, tr = 1, 2, 3, and otherwise. 
Orthogonal bases of £a with respect to ( , )o-.(t for a = 1.2. 3, respectively, are said to be BO bases 

2.3.2. Second generalization. Let £ be an TV-dimensional vector space endowed with two distinct Euclidean 
inner products ( , ) and ( , )'. For a given Ni- dimensional subspace £\, does their exist a subset T of vectors 
orthogonal to £1 with respect to {,) and (, )' at the same time? It follows from section 2.1 that T is the 
subspace Pi £^ , intersection of both orthogonal complements of £1 , £^ with respect to ( , ) and £^ with 
respect to ( , )'. Let {e'^, . . . ,ejv-jvi} be any basis of £^' . There exist uniquely defined vectors Uj G £\ and 
Sj e £1 such that e'j = uj + £j, j = 1, . . . ,N — Ni. Then, {ei, . . . , sn-Ni} is a basis of £1- and for any 
vector V in J^, one has 

N-Ni N-^Ni N-Ni N-Ni 

"j^J= Y1 ^'i^'i ^ K-"j)£j= X) a'^Uj=0. (2.33) 

i=i i=i i=i i=i 

The last equation above vanishes since the left-hand side is a vector in £^ while the right-hand side is in 
£\. This implies that the dimension of is N — N\ — r, where r is the rank of {u\, . . . , un-n^}- (This rank 
is intrinsic, i.e. it only depends on £^ and £1 , independently of the choices of bases for these subspaces.) 
These N — Ni u- vectors belonging to the iVi-dimensional subspace £\, one has r < inf(A^i, — Ni). As 
a first example, if these it- vectors are linearly independent, r = N — Nx (requiring N — Ni < Ni) and the 
dimension of vanishes, i.e. J-" = {0}. As a second example, if ff*" = ff*" , although the two inner products 
are distinct, all the m- vectors vanish, r = and !F — £^. 

Remark. Most of the properties considered in this section 2 can be extended to Hcrmitian vector space 
over the complex field and also to infinite-dimensional Hilbert spaces, the typical spaces used in quantum 
theory 



3. Application to vector spaces of polynomials 

This study can be extended, with similar conclusions, to p > 2 subspaces £0^ and inner products ( , )a.p+i, 
not all the same, defined on £ for cr = 1, . . . ,p. Then, the question is to determine the subspace £p+i such 
that £ = £cr and {£a , £p+i)cr,p+i =0, a = l,...,p. 

See, e.g., [8] chapter V, section 2 and chapter VII, section 4, [11] chapter 3 or [19] chapter VI. 
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3.1. Euclidean vector space of polynomials 

Let be a non-constant and non-decreasing real function on the real domain V such that all the moments 
are finite, 



dij{x) x"- < oo n = 0, 1,.... (3.1) 



If /i is absolutely continuous d^{x) = w{x)dx where, w is a nonzero and non- negative weight function. The 
associated inner product of real functions / and g in the class of square integrable functions L'^^iV) is denoted 

by 

{f,g):= [ dnix) f{x)g{x) . (3.2) 
Jv 

This inner product is symmetric, {f , g) = {g , f), and it follows from the hypotheses on jj that it is positive 
definite, i.e. (/,/)> except for / = 0, where it vanishes. Now, for N a positive integer, let Vn be the 
iV-dimensional vector space over the real fields of real polynomials of degree at most A'^ — 1 (by definition, 
Vo := {0}) Any inner product defined as above endows Vn with an Euclidean structure. Then, all the 
general results considered in section 2 can be applied. Furthermore, new properties are available in connection 
with (i) the additional characteristics of the vectors, i.e. the degree, possibly the parity or the monomial 
character of polynomials and basically, a richer algebraic structure, the multiplication endowing the set of 
polynomials with the structures of a ring and of an algebra over the real field and (ii) the definition of the 
inner product in terms of an integral. The theory of orthogonal polynomials has been extensively studied 
The G-SO of the monomial basis {a;°, x^,. . . , x'^~^} of Vn provides the unique standard orthogonal basis of 
polynomials {Qo, ■ ■, Qn-i} such that Qn is a polynomial of degree n with a given arbitrary real nonzero 
factor kn as coefRcient of a;" , 

{Qj , Qk) = hj6j,k and Qj = kjX^ + 0{x^-'^) . (3.3) 

To avoid the one unit lag between the index of these basis vectors and the degree, from now on, the indices 
indicating the degree will start at instead of 1. This induces some trivial changes in the general formulae 
(2.2) (2.10) which are given in appendix A. 

Classical orthogonal polynomials (i.e., Hermite, Laguerre and Jacobi polynomials) correspond to particu- 
lar measures {V, dfj.}. These polynomials arise frequently and have been studied in great detail Particular 
SBO polynomials associated with these polynomials (as those considered in [1, 3] for the Hermite case, and 
in [2] for the Laguerre case) will be studied in forthcoming papers. 

3.2. Standard block orthogonal polynomials for two subspaces 

Let us apply the general considerations of sections 2.1 and 2.2 to the following case: (i) £ is the vector space 
of polynomials Vn with three Euclidean inner products (,),(, )i and ( , )2 defined as in equation (3.2) by 
the measures {V,d^}, {Vi^d^i} and {V2,dfi2}, respectively; (ii) £i is any subspace oiVw, e.g., defined by 
a,ny sequence {ei, . . . , ejVj} of A'^i linearly independent polynomials in Vn- For £i = Vq or Vn, the problem is 
trivial. For < A'^i < iV, the general procedure given in section 2.2 can be applied to get BO bases of £i and 
£2 satisfying the orthogonality relations (2.16). The first step is to determine the unique subspace £2 = £1 
complementary and orthogonal to £i with respect to the inner product (, ). For that purpose, the basis 

Let us recall the basic property: the Vj, j = 0, 1, . . . are nested, i.e. < j < fc, Vj C Vk- 

See, e.g., [9] section 35 or [10] chapter V for a formal study of polynomials. 

The standard textbook on this subject is [12]. See also, e.g., [13] and [14] chapter X. 

There is an infinity of orthogonal polynomials, see footnote 22. As an example, G-SO of the monomial 
basis taken in the reverse order, i.e. {e\ := x^~^ , . . . ,eN '■= x°}, generates the unique orthogonal basis 
of polynomials {Ri := Ei,. . . ,Rn := En} such that Rj is a polynomial of degree at most N — 1 and 
in which the lowest degree monomial is x^^^ with a given arbitrary nonzero aN-j,j as coefRcient: Rj = 
aijx^~^ + a2,jX^~^ H h aN-j,jX^~^ , thereby defining another kind of standard orthogonal polynomials. 



^2 See, e.g., [14] sections 10.6-10.13, [15] chapter 22 or [16] section 8.9. 
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{ei, . . . , cni} of £i is completed to get a basis {ci, . . . , ejv} of Vn- Tlie G-SO of this basis provides a basis 
{El, . . . , -Ejv} of orthogonal polynomials with respect to the inner product ( , ). Then, £^ is defined uniquely 
by the sub-basis {Ejy^+i, . . . , E^}. Finally, G-SO of {ei, . . . , bn^} (or of {Ei, . . . , -Eati}) with respect to ( , )i 
and of {En-^^i, . . . , E^} with respect to ( , )2 provides BO bases of £i and respectively. It should be 
noted that this general procedure applies whatever the degree of the polynomials defining £i is. 

Henceforth, let us focus on the special case where, £\ = Vi, i.e. the subspace of polynomials of degree 
less than i, for given < i < N . One basis of this subspace is {x^, . . . ,x'^~^}. Then, following the procedure 
above, the unique subspace V^i^ := £i~ such that, 

{'Pi , Hn) = and VN=Vi® V^j^ (3.4) 

is defined uniquely by the sub-basis {Qi, • • • , Qn~i} of standard orthogonal polynomials with respect to the 
inner product ( , ). For given arbitrary real nonzero finite fc„ := b^\, n = 0, . . . , N — 1, these polynomials, 
satisfying equation (3.3), are defined uniquely by equations (A.l)-(A.IO). 

The inner products ( , )i and ( , )2 induce a new Euclidean structure on Vn with the inner product ( , )o 
defined as in equations (2.13) and (2.14), 

w , ^ ^ {P,Q)0-= {Pi;l , «i;l)l + {Pi;2 , qi;2)2 ■ (3.5) 

The BO bases {Pifi, . . . , Pi-i-i} of Vi and {Pi,i, . . . , Pi-N-i} of V^j^i form an orthogonal basis of Vn with 
respect to the inner product ( , )o such that, as in equation (2.16), 

{{Pi-.m 7 Pi\n)l — E^i^m^m^n ^7 ^ — 0, . . . , i 1 
{P^■,m , Pi;n) =0 m = 0, . . . , i - 1 n = i, . . . , N - 1 (3.6) 

(,Pi;m 5 Pi\n^2 — Hi'^^Sm,n ^? ^ — ^5 • * • ? ^ 1 • 

Riquiring that Pi-n be a polynomial of exact degree n, the equations above determine uniquely these polyno- 
mials up to an arbitrary nonzero factor in each polynomial: 

- the Pi;„, n = 0, . . . ,i — 1, which only depend on the inner product (, )i, are the standard orthogonal 
polynomials for the measure {X'i,d/tii} as given in appendix A, with nothing new. In what follows we will 
not be concerned with them; 

- the Pi-n, n = i,. ■ ■ ,N — 1, which depend on i (regardless to the inner product ( , )i) and on both inner 

products ( , ) and ( , )2, arc the standard block orthogonal (SBO) polynomials. They can be obtained using G- 
SO to orthogonalyse the basis {Qi, . . . ,Qn-i} with respect to the inner product ( , )2. Indeed, the procedure 
preserves the degree of the polynomials and one has 

n n 

Pi;n — ^ ^ Qm-'^i;m,n Qn ~ ^ ^ Pi;m-^i;m,n ^i\n,n ~ ^i-n,n ' (^*'^) 
m—i m=i 

Hence, with Pi-n a monic polynomial, 

Pi',n • — ^i;n Pi\n Pi\Ti • — ^ ~t~ P4',n-^ ~t~ ^i\n^ ~t~ 0(iC ) (3.8) 

and equations (3.7) and (A. 3) yield 

^i',n ~ ^n-^i;n,n ~ {bn^n-^i\n^n) ('^*^) 

Ri;n = r-„ -h — (3.10) 

-^i;n,n 

^ , kji—l -^iin—l.Ti ^ , kn—2 ^z;n— 2,n /o i i \ 

Si-,n = Sn + -r -2 Vn-l + -7 T. " (^-H) 
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Then, given arbitrary nonzero finite Bi-jj (or equivalently Ki-j) for j = . . . , TV — 1, and with 

7j,fc := {Qj , Qk)2 J, = . . . , iV - 1 
it follows from equations (2.19) and (2.4)-(2.10), that for i<m<n<N — 1^ 

1 



p. 



■det 



(7i,fc) 

(Qk) 



1 j,k=i,...,n 



Zi-i-i := 1 



TT. ff. jL, 

K A- 



1 



det 



{Sm,k) 

(7i,fc) 
{ln,k) 



(3.12) 

(3.13) 

(3.14) 
(3.15) 

(3.16) 
(3.17) 



For j = 0, the procedure above amounts to orthogonalize the basis {Qq, . . . , Qn-i} of Vn with respect 
to ( , )2- Then, Po;n is the unique standard monic orthogonal polynomial Q2;n of exact degree n, defined for 
the measure {V2,dfj.2}- Prom equation (A. 6), one has 



det 



[(Qk) 



det 



J— 0,...,n — 1 
k=i,...,n 



Zo-n- 



Q2;n : = 



(C2;j+fe) 



j = 0,...,: 
k = 0,.. 



det [c2;j-+fc] 



j,fe— 0,...,n— 1 



where 



dii2{x) a;" 



(3.18) 



(3.19) 
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Thus, for i = 0, the general results are still valid with Vo := {0} and 'Pq^jv •= ^Jv- For i = N, the problem 

is trivial with V^.pf = {0}. 

Finally, in the i,n-plane, the SBO polynomial P^n is associated with the integer coordinate point i,n 
such that <i <n. These points are located in the sector between the two boundary half straight lines: the 
diagonal z = n > 0, along which Pn.n = Qn] and the y-axis i = 0, n > 0, along which Po;„ = Q2;n, where, 
Qn and Q2-n arc the standard monic orthogonal polynomials of exact degree n for the measures {V, dfj,} and 
{'D2,diJ2}, respectively. At the origin, Pq-q = 1. 

Prom equations (3.7), (A.3) and (A. 4), an expansion of Pi-n in terms of the monomials a;™, m = 0, . . . , n, 
as well as the inverse expansion, can be written. As an example, inverting the finite sums over m and £, one 
has 

n m 
m=i e=0 



Pi;n — ^ ] 0'i,mj-^i;m,n — Ci-^^^n 



^=0 



E 

m=sup(i,^ 



(3.20) 



However, we have not been able to get a close form for these connection coefficients Ci-e^n, even in the special 
cases associated with the classical polynomials considered in forthcoming papers. 

3.3. Properties of standard block orthogonal polynomials for two subspaces 

3.3.1. Projection operators. Using the Dirac notation, the general formulae given in sections 2.1 and 2.2.2 
apply readily with only some slight modifications in the notations. In the space Vn^ the orthogonal bases 
of the inner products ( , ), ( , )2 and ( , )o are {Qo, • • ■ , Qn-i], {02;0, • ■ • , Q2-N-i} and {P.-q, Pi-N-i}, 
respectively. Prom equation (3.18), Q2;n oc Po;n, therefore, {Q2;0, • • • , '32;jv-i} can be replaced by the basis 
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{Po;o, • • ■ ,Po;JV-i}- Thus, the orthogonality and the closure relations for these orthogonal bases read for 

0<'i<iV, 

{Qj\Qk) = hj5j,k 2{P0;j\P0;k) = Ho-jSj,k o{Pi-j\Pi;k) = Hi,j5j,k j, k = 0, . . . , N - 1 (3.21) 

N-1 N-1 N-1 

E IQ^) V' (Q^l = E H^.} 2{Po;,\ = E 1^^;^) ^i^J " I = " (3-22) 

As chosen in section 2.2.2, the basis {Po;0; • • • ) -Po;iv-i} is obtained from {Qo, ■ ■ ■ , Qn-i} using the G-SO 

procedure with respect to (, )2. (What is peculiar here, is that all the basis vectors Qn, Po-.n and Pi-n arc 
polynomials of exact degree n.) Therefore, with equation (3.7), one has using the closure relations within 

"Pn+l, 

n 

\Po;n) = IVn+i \Po;n) = E {Qm\Po;n) {Qm\Po;n) = ^0;m,n (3.23) 

m=0 
n 

\Qn) = IVn+i \Qn) = E 1-^0;™) -^0;m 2{Po;m\Qn) -f^O;m 2{Po;m\Qn) = Bo;m,n (3.24) 

m=0 

corresponding to equations (2.25) and (2.26), respectively. In other words, {Qm\Po;n) and 2{Po;m\Qn) vanish 
for m > n, or else, the matrices which relate these different bases are triangular. Following equation (2.27), 
the projectors, lip. onto Vi and IIpj_^ onto V^j^, orthogonal with respect to ( , ), are, 

i-1 i-1 

= E (Qnl = E 1^*;") (3.25) 

n N-1 i-1 

= E E (E Bo;j,nAo;n,k)H^l 2{Po;k\ (3.26) 

j=0 k=n n— 

N-1 N-1 

= E = E 1^';") 0{Pi;n\ (3.27) 
n=i n=i 
N-1 

= Iv, - Hp, = E l^0;„) ifoTn 2 (Po;„ | - Hp, . (3.28) 
n=0 

As an example, for i = 1 one has 

Af-l 

n-p, = lOo) /lo ' (<5o| = |fo;o) -Bo;o,o E Mo,kH^l 2{P0;k\ • (3.29) 

fe=0 

Let lis emphasize again that all these equations have to be handled with care since different inner 
products arc involved. Moreover, when the action of the bras ( | and 2( | are expressed with a single integral. 

This kind of relations are given in [2] (9)-(ll) in the special case i = 1. Unfortunately, these equations 
are wrong since the different inner products involved have not been taken properly into account. This is 
especially the case when writing = X]n=i \wn){wn\ which should correspond to our expression (3.25) 
in terms of o( |, not easy to handle as underlined below. Furthermore, (r|fTjv) defined with (9) should be 
exp(— 1/r) times a zero degree polynomial, i.e. a constant corresponding to our \Qo) (in addition, (zn) for 
the Laguerre polynomials should be (—1)" 2 instead of 2). Thus, |crjv) corresponds to some polynomial \Rn) 
of degree N. Nevertheless, it happens that in (10) the bra {aN\ do corresponds to oc (Qo| and moreover, the 
normalization coefficient is such that {Qo\Rn) = 1, yielding (11). Yet, although \Rn){'7n\ is a projector, it 
is not an orthogonal projector. 
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e.g. y f,p G Vn {.f\p) '■= J-p d4i{x) f{x)p(x), this is no longer true for the bra o(/|. Actually, the only way 
to compute o{f\p) follows from the definition of ( , )o by equation (3.5), i.e. it is to split both / and p into 
their components in Vi and V^^^- Then, this inner product reads as a sum of two integrals. 

3.3.2. Integral representations. The relations between the product of differences, a Vandermonde determinant 
and a polynomial alternant of any polynomials p„ = hnX" + 0(a;"~"^) are well-known With y„ := 
{yo, . . . ,?/„_i}, Ai(yi) := 1 and for n = 2, 3 . . ., one has 



(3.30) 



(3.31) 

z 5 • * • ? Vn • 



l<j<k<n- 

(i) Integral representation of Zi-n- In the integral, 

/:=/ dii2{yi)--- di^2{yn) det[Qk{yj)]^. 

J-D2 JV2 

the product of measures is symmetric and each determinant is anti-symmetric in the variables y. 

Therefore, one of the two determinants can be replaced by one of the (n + 1 — i)! terms of its expansion, say 
YYj=i Qjiyj)-! provided one multiplies the result by {n + 1 — i)\. Thereby, integrating independently over the 
variables yi,. . . ,yn {Vj occurs only in the row j), and using equation (3.12), one gets 

/ = (n + 1 - i)! det[7,-,fc]j,fe=i,...,„ . (3.32) 

With equation (3.14), this yields the integral representation of Zi-n for n = i,i + 1, . . . 

1 2 



Zi- 



■D2 



dn2{yn) '^^^[Qk{yj)]j^k=i,...,v 



For i = 0, replacing Qk by y'^, the same steps as above yield the known integral representation 

det[c2;j+fc] .^^0 n = l—rT\\l dH2{yo)--- ( dlJ,2{yn) ^n+\{yn+\f 



(3.33) 



(3.34) 



where, C2-j is defined by equation (3.19). Therefore, using equations (3.3) and (3.30), one recovers the 
relation between the determinants of the metric tensor matrices for the measure {'D2, diJ.2} in the two bases 
{Qj, i = 0, . . . ,n} and {x\ j = 0,...,n}, 

(3.35) 



Zo-n = det[7,-,fc]^.^^o_ „ = [Yl %J det[c2;,+fc]^. „ . 
(ii) Integral representation of Pi-n- Prom equations (3.3) and (3.30), one has 



n— 1 ^ 

A„+i(y„,a;) = A„(y„) JJ(a; - yj) = — — det 
i=o llj=o«'j 



(QkiVj)) 
{Qk{x)) 



(3.36) 



fc = 0,- 



Now, for j = 0, . . . , n — 1, multiplying the jth row of the determinant above by Qjiyj) and integrating over 
yj the elements of this row with the measure {V, d/i} for j = 0, . . . ,i — 1 and the measure {2?2, ^^2} for 
j = i, . . . ,n — 1, one finds with equations (3.3) and (3.12), 

/ dniyo)--- dfi{yi-i) dii2{yi)--- dfi2{yn-i) ^n+i{yn,x)'[[Qj{yj) 
Jt> Jt> Jt>2 Jv2 j—o 



det 



(7j,fc).=.^...,^,-i (7j,fc)^-v-.- 



,n ■ 



(3.37) 



34 See, e.g., [16] 14.311, [18] section 7.1 or [20] equation (C.3). 

3^ This result can also be obtained from the Gram equality. See, e.g., [18] sections 3.7 and 7.4. 

36 



See, e.g., [12] (2.2.7) and (2.2.11). 
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The Laplace expansion ^''of this determinant according to its i first rows, and equation (3.13) yield the 
integral representation: 



11^=0 % 



Zi;n-1 lYjJo ^3 



n-1 



X / dfi{yo)--- dn{yi-i) I dii2{yi) ■ ■ ■ i rf/i2(?/n-i) A„+i(y„, a;) TT . 
Jv Jv Jt>2 Jv2 

For i = 0, this integral representation reads 

h-A^) = ^ / dfj,2{yQ) ■■■ I diJ.2{yn-i) A„+i(y„, a;) TT Qj{yj) . 

^0;n-l J-D2 JV2 j^q 



(3.38) 



(3.39) 



Now, with the same trick as in (i) above, the measure being symmetric and A„-(_i(y„, x) being anti-symmetric 



in 2/0) • • • ) yn-ii one can replace 11^=0^ QjiVj) by the anti-symmetric function det[Qk{yj)]j,k=o,...,n-i provided 
one divides the result by n\. Then, from equations (3.30), (3.35) and (3.36), one recovers the known integral 
representation for standard monic orthogonal polynomials of degree n for the measure {X>2, dfi2} 



Po;n{x) 



i\ det[c2; 



J+*=J 7,fe=i 



-/ dll2{yo)--- I dll2{yn-l) ^n{ynfW{x ~ yj) 

_1 ^X>2 JX>2 j-=0 



(3.40) 



3.3.3. Domains and measures symmetric with respect to the origin. Let the domains V and D2, and 
the measures d/i and d^2 be symmetric with respect to the origin, e.g., D = [a,b] with a = —h and 
dii{x) = w(x)dx with an even weight function w{—x) = w{x), and similar relations for the measure {T>2.i dfj,2}- 
Then, as the standard orthogonal polynomials Q„, n = 0, 1, . . . (see appendix A), the polynomials Pj;„, n = 
+ . . are even or odd according as n is even or odd, 

fi;„(-ar) = (-l)"Pi;„(a;). (3.41) 

Indeed, Pi;„(— .t), n = i,i + l, . . . satisfies the same orthogonality relations (3.6) as _Pi;„(.x), and since there are 
unique polynomials satisfying these orthogonality relation and the normalization condition (3.8), comparing 
the coefficient of x" yields equation (3.41). Actually, the even and odd polynomials separate from each other 
and equations (3.13)-(3.15) take the form for < i < n 

' (72j,2fe) ■ 



-P2i-l;2n — P2i;2n 
.(e) 



L 



■det 



^i:n-l 



zn := 



det [72j,2fc I 



{Q2k) 

j,k=iy....,n 



4^1 := 1 



H. 



2i-l:2n 



2i:2n 



-2 i\n 



2n (c) 
^j;n-l 



2i;2n+l 



2i+l;2n+l 



2n+l (o) 



■ det 



{l2j + l,2k+l] 
{Q2k+l) 



det[-r2j+i,2k+i]j^k^i^^ 



2i;2n+l 



= H. 



2i+l;2n+l 



= k. 



7(0) 



(3.42) 

(3.43) 
(3.44) 

(3.45) 
(3.46) 
(3.47) 



3'^ See, e.g., [18] section 2.5.1. 

See, e.g., [12] (2.2.7) and (2.2.10) where, the polynomials are normalized according to equation (A. 11), 
or [14] 10.3 (3)-(5). 

^® It is recalled that here, and in what follows, if i = 0, the expressions with the index 2i — 1 do not occur. 



17 



As for the standard orthogonal polynomials (see equations (A. 17) (A. 28)). these relations ean be obtained in 
two ways: (1) it can be checked, as for equation (2.6), that these monic polynomials satisfy the orthogonality 
conditions (3.6) for ( , ) and ( , )2 which define them uniquely; (2) taking advantage of the checkerboard 
structure of the matrix associated with the metric tensor, since 7^ ^ vanishes ii j + k odd, the determinants 
which occur in the general formulae (3.13)-(3.17) can be evaluated from the lemma given in appendix A (see 
equations (A.13)-(A.15)). As an example, let us compute i-2i;2n, one has 



P2i;2n — 



1 



'2n 



- 

"'Or, 



Z2i;2n-1 
1 



det 



2n V 
^2 



■det 



(7j,/c) 

L (Qk) 

(72j,2fc) 
{Q2k) 



j = 2i,...,2n 
k = 2i,...,2n 



det[72j+ 



1.2fc+ljj-fe=i,...,„_l 



(3.48) 



^2i;2n-l 

which, as expected, does not depend on the odd polynomials Q2j+i, j = i, ■ ■ ■ ,n — 1, and .^2i;2n-i reads 



Z2 



i:2n-l 



= det[7j^ 



j,fc=2i,...,2n- 



-1 = det [72j,2fc]^.fe^. det[72j+i,2fe+i]^. j^^. 



,,n-l ' 



(3.49) 



Thereby, one gets the expressions (3.42) and (3.43) of -P2i;2n and Z^?^, respectively. 
It follows from equations (3.42) and (3.45) that 

if (-1)^+" = 1 then Pi_i;„ = Pi;„ < i - 1 < n . 



(3.50) 



In other words, in the i,n-plane, two neighboring points having the same ordinate n are associated with 
equal SBO polynomials if both coordinates i and n of the right-hand point have the same parity. As a special 
case, the points immediately on the left of the boundary i = n > 1 are also associated with known SBO 
polynomials, since from equation (3.13) one has 



— l;i — Pi;i — Qi ^ ^ 1 . 

Therefore, in addition to Pq-q — 1, one has Pq-i — Pi-i = x. 

Note that if i = 1, setting i = in equation (3.45) yields -Pi;2n+i for n = 0, 1, . . ., 

-Pl;2n+1 = Q2;2n+1 = Po;2n+l 



(3.51) 



(3.52) 



as it should be from equations (3.18) and (3.50). Actually, the identity between Pi-2n+i and Q2:2n+i follows 
readily from parity arguments. Indeed, the odd standard orthogonal polynomials with respect to ( , )2, i.e. 
Q2;2n+i) n = 0, 1, . . ., are orthogonal with any even polynomial with respect to any inner product defined 
with an even measure. Thus, they are orthogonal to with respect to ( , ) and orthogonal to Pi-2n with 

respect to ( , )2. 

3.3.4- Normalization. As for the standard orthogonal polynomials (see appendix A) three particular choices 
of normalization can be considered with interest: 

(i) Ki-^n = 1, then Pi-n = Pi;n is a monic polynomial; 

(ii) Hi-n = 1, then the polynomials Pi-n, n = 0, 1, . . . are orthonormal. In addition, choosing K^„/kn > 0, 
from equation (3.15), one has 



— = ( ^ ) Pi;n = {Zi-n-lZi-n) ^ det 



{Qk) 



p. . — 7~^l'^n. ■ 

^ i\i — ^iA 1 



(3.53) 



(iii) Ki-^n = knZi-n-i, then Z-factors in denominators remain only in equation (3.17), in particular. 



^i;n — Z-i-n — lZi-^^ Pi;n — det 



(7j,fe) 

(Qk) 



(3.54) 



j = i,...,n — l 
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3.3.5. Linear recurrence formula with respect to the degree n. Standard orthogonal polynomials fulfill the 
three term linear recurrence formula for n = 0, 1, . . . 

Qn+\ = {AnX + Bn)Qn - CnQn-\ ^" = -7-^ -Bn = - r„) ^n = ^-^-^^- (3.55) 

Following the same steps, with Bi-^ ■= Ri-n+i — Ri;n, the polynomial Pi;n+i — {x + Bi-^n) Pi;n is of degree at 
most n — 1. Therefore, it reads X]j=o '^^■3 where, from the orthogonality condition (3.6), 

Oii-jHi-j = {Pi-j , Pi;„+1 - (a; + Bi.ri)Pi;n)Q = - {Pi;j , xPi-n)^ ■ (3.56) 

But now, for i < n, (Pi-j , a;Pj;„)o is not equal to {xPi-j , Pj;n)o, since the inner product ( , )o is not defined by 

a single integral. From equation (3.5), to evaluate {Pi-j , xP^n)o, one has to split xPi-^i into its components 
in Vi and in V^^, e.g., by expanding xPi-n on the basis polynomials Pi-j, j = 0, . . . , n + 1. This can be done 
as follows using equations (3.7) and (3.55), for i < n, 

n n 

^Pi;n ~ ^ ^ ^Qm-^i-jm^n ~ ^ ^ -^m {.Qm+1 ^mQm ~t~ ^TrtQm—l^'Ai;m,n 
m=i m=i 

n+1 n+1 

= Pi;i-lQi-l + ^ Pi;mQm = Pi;i-lQi-l + ^ 'ni;mPi;m (3.57) 



where, f^i-.i-i, Pi-m and rji-^m for m = i, . . . , n + 1 can be easily written in terms of Ai-j^k, Bi-jj^, Aj, Bj and Cj 
ioi i < j < k < n + I. In the subspace Vi, the expansion of Qi-i in terms of Pi-j, j = 0, — 1 depends on 
the inner product ( , )i In the generic case, all the coefficients ai.j = —rji-j, j = i, . . . ,n — 1 are nonzero, 
while for the standard polynomials, only the equivalent coefficient for j = n — 1 docs not vanish. Thus, 
there is no generalization of the standard three term linear recurrence formula involving a constant number 
of terms. 

However, five term linear recurrence formulae arc derived in a forthcoming paper II for special cases 
associated with Hermite and Laguerre polynomials. These relations, as also some differentiation formulae, 
are based on differential equations satisfied by the weight functions w and W2 we consider. 

3.3.6. Linear recurrence formula with respect to both i and n. Equations (3.7), (3.16) and (3.17) allow to 
relate polynomials P,;„ corresponding to different values of i. With < i < j < n, inverting the finite sums 
over m and £, one has 

n m n n 

Pj-,n='^^(^^Pi\eBr,e,rnjAj.m,n=^^Pr,eCiJ-e,n Cij-e^n '■= Bi-i^rnAj-m.n (3.58) 

m=j t=i e=i m=sup{j,i) 

where Ci^i-e,n = Se,n- Furthermore, since Pt-n and Pj-n are monic polynomials, Cij.n,n = 1- Setting j = i + 1, 
and provided one is able to compute the coefficients Ct^i+i-e^n, the linear recurrence formula (3.58) allows to 

determine step by step Pi-^ for < i < n, starting with Po:„ given by equation (3.18). This is done explicitly 
in a forthcoming paper II for special cases associated with the Hermite and Laguerre polynomials. 

3.3.7. Properties of the zeros. Several properties of the zeros of standard orthogonal polynomials have been 
known for a long time: e.g., all the zeros are real, simple and located in the support T> := [a, b] of the measure 

*° See, e.g., [12] theorem 3.2.1 or [14] 10.3 (7) and (8). A„ and P„ are such that Qn+i - {AnX + P„)<3„ 
is of degree at most n — 1, thus, it reads J2]=o '^jQj where, from the orthogonality relation (A. 2), for 
j — 0, ... ,n — 1, ctjhj = —{Qj , AnxQn) = —An{xQj , Qn)- Since xQj is of degree j + 1, the last inner 
product vanishes except for j = n—1. The coefficient of a;" in xQn-i being kn-i, this fixes a„_i = — C„ as 
given in equation (3.55). 

If ( , )i = ( , ), then Pi-j = Qj, j = 0, . . . ,i — 1 and a^j vanishes for j = 0, . . . , i — 2. 
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and also, the zeros of polynomials with consecutive degrees separate each other Wc saw previously that 
the SBO polynomial Pi-n coincides with standard orthogonal polynomials in several special cases, i.e. from 
equations (3.13) and (3.18), 

Pn;n = Qn Po;n = Q2;n n = 0, 1, . . . (3.59) 

and furthermore, if the domains V and V2, and the measures rf/x and dfi2 are symmetric with respect to the 
origin, from equations (3.51) and (3.50), 

Pn-l;n = Qn Pl;2n-1 = Q2;2n-1 U = 1,2, . . . . (3.60) 

Now, what about the zeros of P^n in the remaining generic case < « < n? Starting the proof as for 
the standard orthogonal polynomials (see footnote 42), one has from equation (3.6), 

{Pi;n , Pi;0)0 = {Pi;n , Pi;o) OC / dll{x)Pi,n{x) = . (3.61) 

Hence, the measure being positive, there exists at least one value of x in D := [a, b] where Pi-n changes sign. 
If Pi-n changed its sign in X> at m < i points Xi,. . . ,Xm, the polynomial Rm ■= YljLi{x — Xj) of degree m 
would be such that, 

iPi;n , Rm)0 = {Pi;n , Rm) = I dn{x)Pi.n{x)Rm{x) = . (3.62) 

Jv 

This contradicts the fact that Pi.„(x)i?m(x) would have a constant sign throughout V. Therefore: for 
< i < n, Pi-n ^as at least m > i distinct real zeros of odd order in V. Also, m < n, since P^n is of 
exact degree n. Since the coefficients of the polynomials considered are real, if a: is a zero, so is its conjugate 
complex. Therefore, it follows that P„_i;„, n = 1.2,... has n real and simple zeros in T), generalizing 
thereby the property following from equation (3.60) only for symmetric measures. On the other hand, the 
usual argument above cannot be applied when i <m<n since, although {Pi-n , Rm)o = 0, this inner product 
is no longer defined by a single integral. 

Concerning the relative positions of the zeros of SBO polynomials with a given i and consecutive degrees, 
we have not been able to prove a general property. Let us only recall that the proof of the interleaving of 
the zeros for standard orthogonal polynomials is a consequence of the Christoffel-Darboux formula which 
itself follows from the three term recurrence formula. It has already been emphasized in section 3.3.5 that 
such a recurrence formula no longer holds for SBO polynomials. 

4. Conclusion 

In order to take into account the particle-number conservation within the density functional theory, B Giraud 
et al [1-3] considered a new sets of real 'constrained polynomials'. These polynomials satisfy a constraint of 
vanishing average with a positive measure and are orthogonal with respect to a inner product defined by a 
second positive measure. The linear constraint considered can be viewed as an orthogonality relation with a 
constant polynomial, with respect to another inner product defined by the first measure. This allows us to 
recast the determination of these polynomials into a more general problem of finding particular orthogonal 
bases in an Euclidean vector space endowed with distinct Euclidean inner products. 

Recalling basic properties of linear algebra, it is shown that, given an Euclidean vector space £ with 
the inner product (, ), and any i-dimensional subspace Si, there exists a unique complementary and or- 
thogonal subspace of codimension i. Endowing each of these subspaces with different Euclidean inner 
products (, )i and (, )2, respectively, yields to define Si and S;^ as BO subspaces. This induces a new 
Euclidean inner product ( , )o on S, which coincides with (, )i on Si and with (, )2 on S^^, and vanishes 
otherwise. An orthogonal basis with respect to ( , )o is called a BO basis, providing a general frame to 
study particular constrained polynomials. We give a general strategy to determine such a basis using a two 

^2 See, e.g., [12] theorems 3.3.1 and 3.3.2. 

^3 See, e.g., [12] theorem 3.2.2 or [14] 10.3 (10) and (11). 
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step G-SO procedure. Let {eo, . . . .ch-i} be any basis of A first step determines S^^ by G-SO of any 
basis {eo, . . . , Cj-i, Ci, . . .} of with respect to the inner product ( , ). This provides an orthogonal basis 
{£^0) El, . . .} oi £ such that {Ei, E'i+i, . . .} is a basis of In a second step, a G-SO of {Ei, Ei+i, . . .} with 
respect to the inner product ( , )2 provides an orthogonal basis {(pi, 4>i+i, . . .} of £j^. Such a basis depends 
on £i and on the two inner products ( , ) and ( , )2- (The standard problem of getting an orthogonal basis 
{(po, ■ ■ ■ , <Pi-i} of £i with respect to ( , )i brings nothing new.) All the relevant quantities in the G-SO (e.g., 
connection coefficients between bases) are expressed in terms of determinants with metric tensor components 
as entries, except possibly for the last row. Furthermore, the G-SO is a step by step procedure well-suited 
for numerical calculation. The importance of the non-degenerate and/or the definite character of the inner 
products considered is underlined. Since different inner products are considered simultaneously, the Dirac 
notation must be handled with care, in particular to express the orthogonal projectors onto £i and £:[- . All 
these considerations can be extended to Hermitian vector spaces and to Hilbert spaces. 

These general results apply readily to the special case of vector spaces of real polynomials with Euclidean 
inner products defined with positive measures. The generalization of the problem of constrained polynomials 
is as follows. Let (, ) be an Euclidean inner product and {eo, . . . , ei_i} be any i linearly independent 
polynomials defining the i-dimensional subspace £{ . The polynomials P, satisfying the i constraints (ej , P) = 
0, j = 0,...,i — 1, belong to the unique subspace £^ of codimension i which is complementary and orthogonal 
to £i with respect to ( , ). An orthogonal basis of £^ with respect to a distinct Euclidean inner product ( , )2 
defines a BO basis. Now, the structure of polynomial algebra with the multiplication law x'^ = x^'^'^, and 
the definition of the inner products in terms of integrals provide additional properties. 

We make a systematic study of these properties in the special case where £i is the subspace Vi of real 
polynomials of degree less than i. Then, the first step to define T'/" coincides with the determination of 
the standard orthogonal polynomials Qo,Qi,. . . with respect to the inner product ( , ), Qn being of exact 
degree n. Requiring that the BO basis polynomials Pi-n '■= (pn, n = i,i + 1, . . . of V^, be of exact degree n, 
defines them uniquely apart from a constant multiplicative factor for each polynomial. We call them the SBO 
polynomials. Emphasizing the similarities to and the differences from the standard orthogonal polynomials 
(basically due to the fact that inner product ( , )o is no longer defined by a single integral), we investigate 
the projection operators, the integral representations, the special case of symmetric measures with respect 
to the origin, the normalizations, the recurrence formulae with respect to the degree n and with respect to 
both i and n and finally the zeros. For physical applications, it is useful to consider the special case where 
the inner product ( , ) is defined by the classical measure corresponding to a classical orthogonal polynomial 
set. This will be done in a forthcoming paper II for Hermite and Laguerre polynomials. 

Applying the general formulae, a similar calculation can be done with a choice of the subspace of 
polynomials £i different from Vi, Then, the BO polynomials are not what we call the standard ones, and 
a priori, a new similar study starting from the general formulae (2.4)-(2.10) is required. Note that if the 
measures considered are symmetric with respect to the origin (see section 3.3.3), and if £i is a subspace of 
polynomials with a given parity, then, due to trivial parity arguments, the monic BO polynomials of the 
opposite parity coincide with the monic standard BO polynomials Q2;n of the same degree. For example, 
for a given m > 0, let £i be spanned by {.t^"'+^} associated with i = 1 constraint of zero odd (2m -|- l)th 

moment. Then, the even monic BO polynomials read P2n = Q2;2n, n = 0, 1, The odd BO polynomials 

require a new determination, except for m = 0, where, the even and odd polynomials being decoupled and 
X being the lowest degree odd polynomial, the problem reduces to the same study as for the odd SBO 

polynomials in the case i = 2. Thus, one has -Pj^+i = f2;2n+i) n = 1,2, Such a case has already been 

considered in [1] 

More generally, instead of constraints all associated with the same inner product (, ), one may think 

The 'Hermite polynomials constrained by a zero momentum', considered in [1] section 3 (17), with the 
weight functions w := e~^^ and W2 ■= , is not the case i = 2. Indeed, the one-dimensional subspace £i 
considered there, and spanned by {x^}, is not 7^2 spanned by {x°,x^}. Nevertheless, it follows from above 
that the monic BO polynomials are: P2„ = Q2;2nj n = 0,1, . . . and P2n+i = P2;2n+i, n = 1,2,.... However, 
this simplification would no longer hold for measures not symmetric with respect to the origin and/or, e.g., 
Hermite polynomials constrained by zero first and second moments, i.e. now with £2 spanned by {x^,x^}). 
Then, the BO polynomials would not be the standard ones. 
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to use a similar approach to take into aceomits constraints associated with different inner products. As a 
first example of generalization, does there exist vectors v (e.g. polynomials) such that, {ei , v)i = and 
(e2 , v)2 = (i.e. for two given distinct vectors ei and 62, and two given distinct inner products (, )i 
and ( , )2)? We show in section 2.3.1 that this generalization of BO subspaces to more than two subspaces 
may have no solution, a unique solution or an infinite number of solutions according to the inner products 
involved between pairs of subspaces. As a second example of generalization, does there exist vectors v such 
that (ei , w) = and (ei , f )' = (i.e. for one given vector ei, and two given distinct inner products ( , ) 
and ( , )')? We argue in section 2.3.2 that, for this example, v belong to a subspace of dimension 1 or 
according to the inner products considered. One of the referees brought our attention to a similar study 
on multiple orthogonal polynomials or Hermite-Pade polynomials [21], e.g., with < m < n, polynomials 
Pm,n of degree n orthogonal with respect to a weight function wi to , j = 0, . . . ,m and orthogonal with 
respect to another weight function W2 to x^ , j = 0, . . . ,n — m — 2. Assuming the weigt functions are classical, 
i.e. solution of Pearson's differential equation {(piu)' + tl)w — 0, with polynomial coefficients cj) and t/i, these 
polynomials have many new applications and nontrivial generalization of the theory of standard orthogonal 
polynomials. 
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Appendix A. Standard orthogonal polynomials 

As considered in section 3.1, let Vn be the vector space over the real field of real polynomials of degree 
less than N with the Euclidean inner product (, ) defined with the measure {P, d/i}. With respect to the 
monomial basis {a;°, . . . , x^~^}, the matrix of metric tensor components is a persymmetric matrix known as 
the moment or Hankel matrix which entries are the moments of the measure defined by equation (3.1) 



gj,k = {x^ ,x^)= / diJ.{x) x^+^ = Cj+k j,k = 0,...,N -1. 
Jv 



(Al) 



A.l. Relations between standard orthogonal polynomials and monomials 

For < n < N — 1 and given arbitrary nonzero finite kj := bj^ , j = 0, . . . , n, according to equations (2.4)- 
(2.10), the G-SO of {x°, . . . , a;"} defines uniquely the standard orthogonal polynomials Qj, j = 0, . . . ,n, such 
that Qj is of degree j with kj as coefficient of x^ , 

{Qj , Qk) = hjSj^k j,k = 0,...,n {A.2) 

n 

s„a;"-2 + 0(a;"-3) (A3) 



m=0 
n 

'^^ See, e.g., [18] section 7.4. This basic property follows from the polynomial algebra, see footnote 29. 
As noted in [22] corollary 6.3.2, one can consider as well gj^k '■= {sj , x'^), where, Sj is any polynomial 



of exact degree j, in order to make easier the forthcoming calculations. 
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where Qn denotes the monic polynomial ' . One has forO<m<n<-A/'— 1 





= det[cj+k\j 


fe=0,...,n 


Qn 


= ^det 


[Cj+k) 

[ (x'^) \ 


h 


— k'^h 

• — '^n 


h - ^ 

— ry 




'— kn Q'm,n 


^m,n ~ 


Tn 


— Q'n—l^n 


Sji = Clji 




— l.-l'? 

• — ^m,n 





j=0,...,n — 1 
fc=0,...,n 



Zn- 



■det 



(cj+k) 
{Sm,k) 



■det 



(Cj+k) 
(Cn+k) 



j=0,...,n — 1 
k=0,...,n 



3=0,...,m-l 
k=0,...,m 



{A.5) 
{A.6) 

iA.7) 

(A8) 

(^.9) 
(^.10) 



Apart from the usual choice of normalization for the classical orthogonal polynomials three particular 

normalizations are commonly used: 

(i) fc„ = 1, then Qn = Qn is a monic polynomial; 

(ii) hn — I, then the polynomials Qn, n = 0,1,... are orthonormal. In addition, choosing fc„ > 0, from 
equation (A. 7), one has 



Zn-i\y^ 



~\ Zn ) 



Qn — {Zn-lZn) 



-1/2 



det 



(cj+k) 



- 7-1/2 . 



j=0,...,n — 1 
k=0,...,n 



Qo = Zq 



(iii) kn = Zn-i, then Z-factors in denominators remain only in equation (A. 10), in particular, 

(Cj+k) ' 



hn — Zn—lZn Qn — det 



Qo = 1- 



(All) 



(^.12) 



j = 0,...,n-l 
k=0,...,n 



A. 2. Domain and measure symmetric with respect to the origin 

If, for example, V = [a, b] with a = —b and diJ,{x) = w{x)dx with an even weight function w{—x) = w{x), 
the polynomials Qn, n = 0, 1, . . . arc even or odd according as n is even or odd. Indeed, .t) satisfies 

the same orthogonality relations (A. 2) as Qnix), and since there are unique polynomials satisfying these 
orthogonality relations with the normalization (A. 3), comparing the coefficient of a;" yields Qn{—x) = 
(— 1)" (5„(.x). Actually, the even and odd polynomials separate from each other: Cj+fe vanishes for j + k 
odd and the determinants which occur in the general formulae (A.5)-(A.10) have a checkerboard structure, 
except possibly for the last row. These determinants can be evaluated with the following lemma: 

Lemma A. Let A := (Aj,k)j.k=ci n-i be an n x n matrix with the checkerboard structure Aj^k = if j + k 

odd, except possibly for the last row j = n — 1. Then, one has 

det A = det B det C (^.13) 

where, for m = 1, 2, . . ., 

n:=2m B := (S,- fe := ^2,-,2fc),- fc=o,...,m-i ^ ■= i^i^k := ^2i+i,2fc+i) ._fe^o,...,m-i (^-l^) 

n := 2m + 1 B := (5,-,^ := ^2i,2fc),. ,=o,...,m ^ ^= i'^j^k := ^2,+i,2fc+i)^. ^=0,...,™-! • (^-l^) 

^'^ The coefficient f„ is denoted r„ by [14], see section 10.3 p 158, where the coefficient k'^ := kn Vn of x^~^ 
in Qn is also introduced. 

See, e.g., [14] 10.13 (4), 10.12 (2), 10.9 (8) (for A ^ 0) and 10.8 (5) for the choice of kn for Hermite, 
Laguerre, Gegenbauer and Jacobi polynomials, respectively. 
^9 See, e.g., [12] (2.2.6), (2.2.7) and (2.2.15) where, £)„ denotes Z„. 

23 



The relation above is still true for n — 1, setting det C := 1 as a null determinant. 

For n even, the proof is as follows. Permuting rows and columns of A in order to collect the zero 
elements into two blocks, one gets 



det A = det 



/ B (O),,fe=o,...,™-i \ 

(0) j = 0,...,m-2 ( A2i + l,2fe+l ) 3=0,...,m-2 

k=0,...,m-l fe=0,...,m-l 



(^.16) 



\ {An-l,2k) fe=o,...,m-l (^"-1.2fc+l) fe=0,...,m-l / 



The Laplace expansion (sec footnote 37) of this determinant according to its first m rows yields equation 
(A. 14) (the result does not depend on An-i,2k, A: = 0, . . . , m — 1). For n odd, the proof is similar (then, the 
result does not depend on An-i^2k+i, = 0, . . . , m — 1). 
Now, this lemma yields for n = 0, 1, . . ., 



Z^^ := det[c2j+2k_ 



0,...,n 



:= 1 



Q 



2n 



a2m,2nX — 



7(0) 
■'n-X 



■det 



-2fe ) 



Qo = i 



fc = 0,.. 



^2n _ ^ Q2mb2m,2n 



h2n 



0'2m,2n — 



b2m,2n 



1 
1 



■det 



.(e) 



det 



{C2j+2k] 
(Sm.k) 

(C2j+2fe) 
(c2n+2*;) 



j=0,...,n — 1 
fc=0,...,n 



J =0,. . . ,m — 1 
fc=0,...,m 



4°) := det[c2j+2k+2 

Q2n+1 



j,k=0,...,n 

2m+l 



^2r^+l 



Ci2m+l,2n+lX' 

m=0 



— Q2m+lb2m+l,2n+l 
m=0 

7(0) 



:= 1 

^ det 



7(0) 

■'n-l 



(c2j+2fc+2) 
(2.2fe+l^ 



j=0,...,n — 1 
fc=0,...,n 



«2m+l,2ra+l 



&2m+l, 



7(0) 
1 

1 



■det 



2n+l 



det 



(c2i+2fe+2) 
. {Sm,k) 
{c2j+2k+2) 
(c2n+2fc+2) . 



3=0,...,n- 
k=0,...,n 



3=0,. ..,rr 
k=0,...,n 



A.17) 
A.18) 

A19) 

A20) 

A.21) 

A.22) 

^.23) 
A.2A) 

A.2b) 
A.2Q) 
A.27) 
A.28) 



These polynomials being unique, the relations above can also bo proven by checking directly, as for equa- 
tion (2.G), that the polynomials (A. 18) and (A. 24) do satisfy the orthogonality condition (A. 2) with the 
normalization (A. 3). 



Appendix B. Examples of block orthogonal polynomials for three subspaces 
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Let us only illustrate section 2.3 with the following polynomial examples for equation (2.29), A^i = N2 = 1, 
N = 2>, £1 = Vi, £1^2 = and £ = V3, with distinct inner products ( , defined by the positive measures 
{2?p_(j, Wp^„{x)dx} for 1 < p < cr < 3. Then, each of the three one- dimensional vector spaces £1, £2 and £3 
are defined, respectively, by the monic basis polynomials, 

Po = 1 Pi=x + 00,1 P2=x'^ + 01,2 X + 00,2 (B.l) 

where, ao,i, ai,2 and ao,2 have to be determined from the orthogonality conditions, 

(^o,A)i,2 = (Po , -P2)i,3 = (Pi, ^2)2,3 = 0. {B.2) 

The first equation above defines uniquely £2 in all cases. The last two equations may have for £3 a unique 
solution, no solution or an infinite number of solutions depending linearly on one parameter. 

(i) If, for 1 < /9 < (7 < 3, 'Dp^CT is symmetric with respect to the origin and Wp^a is even, then one gets from 
parity arguments, 

00,1 = ai,2 = ao,2 = -^i^ (5.3) 

where, Cp^a;n '■= jj, dxWp^a{x)x'^ . Thus, in that case, £3 exists and it is unique. 

(ii) If, 

Vp,„ := [0,(X3) Wp^,{x) := e-^x^^-'^ Zp^, > (5.4) 
one finds using the Euler integral ^'^ , Cp^a;n = ^{zp,<T + n), ao,i = — ^1,2 and. 



1 2l,3 \ f ao,2\ _ f ^^1,3 (^1,3 + 1) 

Z2,3 + ao,i 22,3 (22,3 + 1 + ao,i) J \ ai,2 / V ^2,3 (-22,3 + 1) (-22,3 + 2 + ao,i) 



(5.5) 



In the generic case, e.g., zi 2 = 1, ^2,3 = 2 and 21^3 = 3, the determinant of the 2x2 matrix above is nonzero 
and there is a unique solution for £3. For, e.g., 2:1^2 = 1, -2^2,3 = 2 and 2:1^3 = 4, the determinant vanishes 
while the rank of the augmented matrix is 2. Then, there is no solution. For the two linear equations (B.3) 
be the same, one must have 

f (2:2,3 + ao,l) {Z2,3 - 21,3) + Z2,3 = 0^1 ^2,3 - ^1,3 + 1 = 

I (-22,3 + aO,l) {Z2,3 - 21,3) + '2Z2,3 - Zl,3 + 1 = \ flQ,! = . 

This case is excluded since 2:1,2 = —ao,i has to be positive. 

(iii) In the same case as above, except for the inner product ( , )i_2 assumed to be defined now by any domain 
2?!, 2 symmetric with respect to the origin and any even weight function ■u;i,2j ao,i vanishes by parity and 
equations (B.3) and (B.4) still hold. Therefore, for 2:2,3 — 2:1,3 + 1 = with 2:1,3 and 2:2,3 positive, there is 
an infinite number of subspaces £3 defined by the basis polynomial P2 = a;^ — 2:1,3 (2:1,3 + 1) + ai,2(a; — 2:1,3) 
depending linearly on one arbitrary parameter ai,2. 
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